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I. INTRODUCTION 



In field theoretic investigations one is often confronted by the rather formidable task to 
evaluate the one-loop effective action in some nontrivial background field. Until recently, it 
has not been possible in four spacetime dimensions to evaluate explicitly this renormalized 
quantity (including its full finite part), unless some very special background is chosen or a 
priori arbitrary parameters (e.g., mass values) are set to zero. In our recent publications 

m 

we made some headway to this old problem by developing an efficient calculational method 
- a combination of analytic and numerical schemes - for the exact computation of fully 
renormalized one-loop effective actions in radially symmetric backgrounds. For example, 
this method was first applied to the accurate determination of QCD single-instanton deter- 
minants for arbitrary quark mass values producing a result that interpolates smoothly 
between the known analytical massless and heavy quark limits. In Ref . ^ we generalized the 
calculational procedure to calculate the one-loop effective action in any radially symmetric 
background, not just an instanton. In the present paper, which is a sequel to Ref.j^, we 
present some explicit examples and results (including the numerical contributions) to estab- 
lish the efficiency and generality of our method. We also examine the validity of often-used 
approximation methods, such as the large mass expansion and the derivative expansion, 
compared to numerically exact calculations. 

In Ref. we derived some relevant formulas needed in the calculation of the scalar one- 
loop effective action (in Euclidean spacetime), assuming SU(2) background gauge fields of 
the form 

(Case 1): A^{yi) = 2r]^,^aXuf{r)—, (1.1) 
(Case 2): A^{x.) = 2{r]^^iUi)x^g{r)Y, (1-2) 



where /i, z/ = 1, 2, 3, 4, r = y |x| = ^x^x^, rifj_^a (or Vfiiyi) the 't Hooft symbols (3|], and Ui 
a unit 3-vector. Case 1 is inherently non-Abelian, while Case 2 has a fixed color direction 
and so is quasi- Abelian. These backgrounds are characterized by the radial profile functions 
/(r) and g{r), respectively. In each case the spectral problem separates into partial waves 
due to the spherical symmetry. 

Our method has been deliberately developed so that it can accommodate numerical input 
for /(r) and g{r), since this situation often arises in quantum field theory applications. But, 



to illustrate the method more clearly, we choose here specific Ansdtze for the radial profile 
functions. In Case 1, of 'non-Abelian type' (11. II) . we choose the radial function /(r) of the 
form 

m = ^2H{r) , H{r) = ^^L^^ , (1.3) 

with free parameters p and a (under the regularity restriction \a\ > 1). [The BPST instanton 
solution corresponds to the choice a = 1]. In Case 2, of the quasi- Abelian type (11.21) . we 
choose the radial function g[r) of the form 

^(r) = B{l-tanh[/5(/Br-^o)]} , (1-4) 

with three free parameters /?, .^o and B (all taken to be positive). In the limit j3 0, 
(11. 2p then approaches the case of uniform field strength B. For finite /?, (11.21) represents a 
spherical bubble type potential with radius and wall thickness ~ ^75- 

In this paper we calculate the renormalized scalar one-loop effective action (including its 
full finite part) in the gauge field background pertaining to the above two types. The effective 
action is computed for arbitrary choices of the parameters characterizing the shape of the 
background field, not relying on the background field being slowly or rapidly varying, or on 
the particle mass being large or small relative to the scales set by the background field. In 
performing this analysis, we have found that significantly greater calculational efficiency and 
precision can be attained by making a systematic use of higher-order quantum mechanical 
WKB-type approximations 5| for the large partial-wave contributions to the effective action. 
This permits the extension of the high partial wave radial- WKB approximation to lower and 
lower partial waves, and results in dramatic numerical improvements. Finally, we compare 
our results to the predictions based on the large mass expansion and the derivative expansion. 
To our knowledge, this kind of genuinely unambiguous comparison in four- dimensional gauge 
theory has not been made before. 

This paper is organized as follows. In Sec. [TTl we give a short outline of our numerically 
exact calculational scheme and also collect, for later use in the paper, relevant formulas from 
the large mass expansion and derivative expansion for the scalar one-loop effective action. 
Our detailed study on the one-loop effective action with non-Abelian-type backgrounds (II. ip 
is then presented in Sec. IIIII This is followed in Sec. IIVI by the corresponding study with 
quasi- Abelian backgrounds (II. 2p : in this case, the one-loop effective action is essentially that 
of scalar QED. In Sec. |V]we conclude with some related discussions and comments. 



II. BRIEF SUMMARY OF OUR CALCULATIONAL SCHEME AND OTHER AP- 
PROXIMATION METHODS 

The calculational method developed in Refs. can be summarized and streamhned as 
follows. For scalar fields in radially symmetric non-Abelian background gauge fields yl^(x), 
it is possible to express the corresponding (Euclidean) one-loop effective action as a sum 
of individual partial-wave contributions, i.e., J2y=o^ji^]^)y with the J-partial-wave term 
[including the appropriate degeneracy factor] given by radial functional determinants 

,2\ \ 



Tj{A;m)=\n' ^ 



m 



det{-dj 



or, equivalently, by the proper time representation 

ds _5 



Tj{A;m) 



dr tr |Aj(r 



)} 



(2.1) 



(2.2) 



Here, m is the scalar mass, — -Dj {—d"]) denotes the quadratic radial differential operator 
relevant to the J-partial-wave in the given background (or with A^{yi) set to zero), and 
Aj(r, r;s), Aj'^'^(r, r; s) represent the coincidence limits of the related proper-time Green 
functions specified by the radial 'heat' equations: 



i^0+ : 



Aj(r,r';s) = 0, (s > 0) 
Aj(r, r'; s) — > 6{r — r'). 



(2.3) 



(The tildes over differential operators and Green's functions imply that we are here consid- 
ering reduced operators/functions after taking out various radial measure factors 2|]). While 
the quantities Tj are finite individually, the partial wave series J2T=o diverges and should 
be renormalized. We thus introduce the intermediate partial wave cutoff Jl (a large, but 
finite value of J) and express the fully renormalized effective action by two separate terms 



ri.en(^; m) = Tj<j^ {A; m) + Tj^j^ {A; m) 



(2.4) 



where 



Tj<j^{A;m) 
rj>ji(A;m) 



J=0 



f:„ Vdet(-93- 



lim 

A^oo 



ds 



m 



-A2. 



J>Jl 



l^(4^^^(f)/^'^*^^^-^-^ 



(2.5) 



(2.6) 



with 

Fj{s) = dr tr (Aj(r, r; s) - Kf\r, r; s)) . (2.7) 

Notice that the renormahzation counter-term is incorporated in the second term 
rj>j^(A;m); /i is the normahzation mass, and F^^ = i[D^,D^] (with = — iA^) 
denotes the background Yang-Mills field strengths. The first term, Tj^j^, may be evalu- 



ated numerically using the Gel'fand Yaglom method 
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for one-dimensional functional 



determinants. As for the second term 



j>Jl, analytically developed expressions valid for 
large enough Jl were used in Refs. 2j; for this, we made essential use of the quantum 
mechanical radial WKB expansion and the Euler-Maclaurin summation formula. A crucial 
observation is that even though the partial-wave sum in (12.51) is quadratically divergent as 
we let Jl ^ oo, this divergence is exactly cancelled by a similar term originating from the 
second term rj>j^ in (12. 6p in the large J^-limit. Thus, for the sum of the two terms, we 
secure a finite result in the large J^-limit - this is the essence of our numerically accu- 
rate scheme for the determination of the renormalized effective action. The "low" partial 
wave contribution is computed numerically, and the "high" partial wave contribution is com- 
puted analytically in the large Jl limit using radial WKB. The ultraviolet cutoff dependence 
appears in the analytic high partial wave computation, and this allows standard renormal- 
ization techniques to be used, leading to the finite renormalized effective action. For details 
on our renormahzation prescription (and on how to go from ours to other prescriptions), see 
Refs. 111,!, 3. 

Clearly, numerical efficiency of our calculational scheme hinges on how quickly the Jl — >■ 
00 limit is attained; that is, on how much we can lower our partial wave cutoff Jl to 
secure a reliable large- limit value for the above sum of the two terms. (Note that if 
we were able to calculate both parts, i.e., rj<j^ and rj>j^ exactly, their sum would be 
independent of the choice of the cutoff value Jl- An interesting example of this situation is 
the effective action for massless quarks in a single BPST instanton background, where our 



t Hooft's result |3| r(m = 0) = a (i 



technique reproduces exactly and analytically 
-i| - I In 2 + I - 2C'(-1) ^ 0.145873....). In we chose the cutoff value Jrto be rather 
large (of the order of 50), and used the Richardson extrapolation method 10|] to reduce 
the numerical round-off error. There, it was sufficient to use the expression for Tj^j^ with 
O(j^) or smaller terms suppressed. But, in the present, more extensive, study, we have 
found that the computation is greatly improved, both in computing time and precision. 



by including [analytically calculated] higher-order WKB terms in the expression for F j> , 
so that it becomes valid up to the 0{jt) accuracy. Tiese higher-order WKB terms can 
be found straightforwardly from the j-expansion ^ [jjj for the radial proper-time Green 
function. With this procedure, we were able to ensure that the renormalized effective action 
we calculate is independent of the (moderately large) J^-value with relative error of order 
10~^, which is comparable to the final total error involved in our numerical computation 
of (12.51) . In the end, we obtain comparable precision with much lower values of Jl, of the 
order of 10 or 20. This improved precision is demonstrated in sections IIIII and IIVI for the 
backgrounds in (II. ip and (II. 2p . 

We also give the results of some popular approximation schemes for the one-loop effective 
action, which are supposed to capture accurate values in certain limits. First, we have 
the large mass expansion j^, [l2| of the scalar one-loop effective action which is obtained 
most easily with the help of the Schwinger-DeWitt proper-time expansion (or heat kernel 
expansion). Here, from the renormalized effective action, it is convenient to separate /i- 
dependent pieces (and use dimensional considerations) to write 



(Case 1) 
(Case 2) 



i(v4;m) 



g In M 



6 



tr F 

tr F 

(47r)2 



F(mp) 

t( 



m 



(2.8) 
(2.9) 



^ ^BJJ [Any \^B) 
where dimensionless constants are not indicated in an explicit manner. The quantity F(mp) 
or F(^), which is independent of /i, can be expressed for large enough m by an asymptotic 
series of the following form 



F(mp) 



c^^x „ ~(n-3)! 



Here a„(x, x) 



B. 
n 



6 

3,4,5, 



In 



n=3 
2 



[m 



2\n-2 



(47r 



-tra„(x, x) 



(2.10) 



d^x 

\rF 

(47r)2 



E 

n=3 



(n-3)! f d^x 



(m 



2^n-2 



(4vr) 



-tra„(x,x). (2.11) 



denote appropriate coefficient functions in the Schwinger- 
DeWitt expansion: explicitly, for the traces of the 03- and a4-terms, we have jo], [l2| 



tr a3(x, x) 
tr a4(x, x) 



15 
1 



-]:trU^F^xFx,F^.-^iD^Fx^)iD^Fx^) , (2.12) 
o L 15 20 

1 r 1 11 

^ tr ——ri^xPx^r^yPyi^ + 

+ ^F^xFx^F^uF,^ + i^F^x{D,,Fxu){D^F,^) + i^F^x{DxF^,){D^F,^) 



2"!^ FkxFx^F^i/Fi^i^ 



FkxF^i/FxkFu^ ~\~ Fi^xFxuF^pFi 



35 



+—{D^DxF^,){DxD^F,^) 



(2.13) 



where D\F^y = [Dx, F^,^] and D^D\Ffj_y = [D^, [-Da, F^u]], etc. Note that in these expressions 
for a3(x, x) and 04 (x, x), we have not assumed that v4^(x) satisfies the classical equations of 
motion. We also comment that while this large mass expansion is rather simple to use, it is 
in fact an asymptotic expansion, and so its regime of useful applicability is restricted to the 
mass being large relative to ^ or V^, respectively. See Sec. Ullland Sec. IIV]for comparisons 
of the large mass expansions (12.101) and (12. lip with our exact numerical answers. 

There is another well-known approximation method to the one-loop effective action, the 
derivative expansion. The leading order of the derivative expansion corresponds to using 
the Euler-Heisenberg constant field result 14, 15|, but substituting the inhomogeneous 



fields for the homogeneous ones used to compute the Euler-Heisenberg effective action. This 
approximation is very simple to implement, and is expected to be a good approximation 
when the spacetime variation in the background gauge field strengths is sufficiently 'slow' 
so that we may regard their derivatives as small terms in the effective action. Subleading 
derivative expansion contributions can also be computed, but we will not consider them 
here. A systematic study of the validity range of this method is still lacking (although the 
Borel summability properties of the derivative expansion have been analyzed in a nontrivial 



soluble inhomogeneous QED example in 



16|). We remark here that if the background gauge 



fields are genuinely non-Abelian (as in our Case 1), the convergence character of the related 
expansion - the so-called covariant derivative expansion 



171, 



18 



Il9l . I20I] - is less certain. In 



this paper we restrict ourselves to applying the derivative expansion with our quasi- Abelian 
backgrounds (II. 2p only. In that case, the leading term in the derivative expansion is given 



by the Euler-Heisenberg formula [l3|, Q, ll5l] 



^/B J A Jo Vsinh(Eis) sinh(E2s) 6 



1 , f m \ r (i^x 



In ^ / - — rrrtrFf,,, (2.14) 
6 \^)J (47r)2 ^ > 

where -^iEi and ±Z-E2 denote four eigenvalues of the 4x4 matrix F = (F^^), with 

1 / , „ _ 1 



where 



El = -yj'-VJ'^-g^ E2 = -\/j'+VJ'^-g^ (2.15) 

^ = ^trF^,F^„ g = ^tie^uXnF^uFx^. (2.16) 

Consult Sec. [IV]to see how the predictions based on this method fare against the accurate 
numerical calculations. 
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III. NON-ABELIAN BACKGROUNDS 



A. Properties of the Background Fields 

We consider here a family of radial background fields described by f 1 1.11) and fll.3p . which 
resemble the single instanton configuration. The parameter a is chosen to be in the range 
I a I > 1, so that A^iiC) is well-behaved at the origin r = 0. When a takes a negative value 
it is convenient to cast the function if(r) in the form 

E{r) = , \ | . (3.1) 

Note that while our configuration with a = 1 is simply the single instanton solution in the 
regular gauge [l|, by choosing a = —1 the single anti-instanton solution in the singular 
gauge is also obtained. In Fig. [T|, we have plotted the shape of the background function 
H{r) for several values of a's. 

For the gauge field (11.11) . the corresponding field strength tensor Ff^^, is 



F^, = y ( + 2{x^r],xa ' X^V^Xa)^ [2H{H - 1) + vH'] ) , (3.2) 



where H' = -^H{r). Alternatively, inserting the expression (11.31) for H{r), we find 

Note that the field strength tensor F^,^ has a singularity at r = when |a| < 1, and 
the second part of (13. 3p . the part proportional to {x^rjyXa — Xyri^\a)x\, vanishes if a = 1 
(i.e. for the single instanton solution). For these instanton-like radial background fields the 
corresponding classical action is readily evaluated: 

]^jdS^tiFl^ = 12^^ j^'^[r\H'f + AH\H -If] 

= 47r^ , , (3.4) 
|a| 

Note that this result does not depend on the scale parameter and has a minimum value 
when a = ±1, i.e., for the single instanton or single anti-instanton solution. 

For the general vector potential considered here, the topological winding number is 



327r 



1 r d 

—e^,xr J dScti F^^Fxr = - dr —{2H^ - = ±1. (3.5) 




H{r 



(a) 




(b) 

FIG. 1: Plots of the radial profile function H{r) appearing in the non-Abelian gauge field (jl.ip . 
as a function of r/p. We have drawn the cases with a = 1,2,3,5 in (a), and the cases with 
a = — 1, —2, —3, —5 in (b). Note that H(r) behaves like a step function if \a\ becomes very large, 
with the step localized at r = p. 



This means that all the configurations corresponding to positive values of a belong to 
the class of winding number 1, and those corresponding to negative values of a to the 
class of winding number —1. Further, notice that self-dual configurations occur when 
d'^ti F^j^ = 5 / dS^ti F^yF^y, which means = 1, corresponding to the BPST instanton 
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or anti-instanton for which a = ±1. 



B. Large Mass Expansion 

Before delving into the exact calculation of the one-loop effective action, we first present 
the result of the large mass expansion. In our background fields we find for the traces 
of leading coefficient functions a3(x, x) and a4(x, x) (see 02.121) and (12.131) ). the following 
explicit results: 

% . f 10 - 13a' - 33a') 

/o pM5(l + r2")6 LV J 

+r2° (5 + 20a + 22a' + 28a^ + 21a^) + (a -a)] , (3.6) 



J d'^x tr 03 = J 



°° dr r 



8a-5 



c/*xtra4 = / -r « 4r^"(l + a)M 28 + 112a + 105a^ + 188a^ + 107a^ 



p^210(l + r2")* 
r^" f-91 - 273a - 21a' + 309a^ + 618a^ + 1164a^ + 694a^ 



+ i (l337- 5250a2- 3399a^ + 11992a^) + (a^ -a) . (3.7) 

Then, after performing the r-integration, the large-mass expansion for the one-loop effective 
action is seen to take the form (we here give the result for the quantity r(mp), introduced 

in (ESD) 

~ , , (1 + tt^), . X 7r(5 - lOa^ + lla^ - 6aS) 1 ,0 x/ 2 

Tlm mp = ~ \ / In mp + ^ ^^^^ e w / ^ + « ' ^ « " 1 

12|a| 1800a° sm(7r/a) [mp)'' 

^7r(-140 + 35a2- 378a4-317a6 + 120a^) 1 

^ 88200a8 sin(27r/a) (mp)^ ^ ' ^ 

Note that, for \a\ = 1 (or 2), taking the limit \a\ 1 (or \a\ 2) in the right hand side of 
(13. 8p should be understood. This large mass expansion result (13.81) will be compared with 
the numerically determined effective action later. 

We make a comment on the length scale parameter p here. The modified effective action 
r(mp) does not depend on the renormalization mass scale p, and so is a function only of 
the dimensionless combination mp. We may then set the size parameter p = 1 during the 
calculation, without loss of generality, and readily restore it in the final result. 
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C. Numerically Accurate Calculation of the Lower Angular Momentum Part 



We now turn to our accurate effective action calculation based on (12.41) . First consider the 
lower angular momentum part rj<j^, given in (12.51) . Note that, in the present backgrounds, 
the partial waves are specified by the quantum numbers J = as described in 

detail in Ref. [2!. We are working with isospin |, so j = / ± |. Using the notation of 2], the 
radial Hamiltonian, representing —D^ in the given partial wave sector, assumes the form 

^'.--^L) = -|^-^|: + ^'. (3-9) 

4/(/ + l) 



with 



j(j + l)-/(/ + l)-^ 



3^^, (3.10) 



^2 



while in the absence of the background field 

. . _|1 _ ?iL , W . (3.n) 

The radial Hamiltonian is independent of the quantum numbers and /s; this introduces 
the degeneracy factor (2j + 1)(2/ + 1) in the partial wave sum below. 

Having identified the relevant quantum numbers, the lower angular momentum part 
can be written as (here, L serves as our partial wave cutoff) 

V 

1/2 

The ratio of two determinants in (I3.12p is determined, according to the Gel'fand-Yaglom 



^ ^ ( det(7^, , + m2) \ 

■,,,JA, ») ^ g(2; + 1) ^ i2, + 1, in [ J^^j; ] , (3.12) 



method 



by the ratio of the asymptotic values of two wave functions as 

det(Hf^^ + m2) R^oo\^f^^{R) ) ' ^' ' 
Here ipij{r) and '?/^p'"^(r) denote the solutions to the radial differential equations 

{nij + m^)^i,j{r) = 0, (3.14) 

(nf'''' + m^)^l'^^{r) = 0, (3.15) 

which have the same small-r behaviors, i.e., 

r ^ : ^ij{r) ~ r\ i^f^^r) ~ rK (3.16) 
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Note that the solution to fl3.15p . which is the modified Bessel function 

(r) = :^?^±^ , (3.17) 

grows exponentially fast at large r, as do the numerical solutions to fl3.14p for the operators 
1-iij + m?. Thus, numerically, it is advantageous to consider the ratio, il)ij{r)/il)f'^'^{r), which 
stays finite for all r. In fact, since we compute the logarithm of the determinant, we can 
directly consider the logarithm of the ratio: 

S,..)W->n^. (3.18) 

which also has a finite value in the large r limit. This function satisfies the differential 
equation 

f^J''^)\{\ + 2jf^y-^ , (3.19) 

dr^ \ dr J \r hi+iynfir) J dr ^ " 

4/(/ + l) 



U{i,j){r) = Vi,j ^2 

4j(j + l)-4/(/ + l)-3 



+ z^TV-z:^.^ (3-20) 



under the initial value boundary conditions 

5(U)(r = 0) = , ^[,.)(r = 0) = . (3.21) 

Noting that the eigenvalues of the total angular momentum j equal / ± | for a given value of 
/, it is convenient to combine the contributions 5'(; ;_,_i-)(r) and S'(;_,„i ;-)(r), which come with 
the same degeneracy factor (2/ + l)(2/ + 2). With this understanding, it is possible to express 
the amplitude (13.121) in the form 

L 

Tj<j^{A;m) = (2/ + 1)(2/ + 2)P(/) , (3.22) 

/=o,i,i,... 

P{1) ^ 5(,,,+ i)(oo) + 5(,+ i,,)(oo) . (3.23) 

Here S(^^i_^ i-j{oo) and S(^i^i^i-^{oo) denote the asymptotic (i.e., r — > oo) limits to the solutions 

of the differential equations in fl3.19p with the potentials 

4:1(1 + 1) H(r) H(r)(H(r)-l) 
n^+i/2(0 = ^^ + (4/ + 3)^ + 3 ^ , (3.24) 



and 



W)^i^ + («.3)(1^.3m_im , ,3.25) 
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respectively. 

Note that the potentials V;^/+i/2(r) and V;_|_i/2,/(r) have the same form except that H{r) 
in one expression gets replaced by (1 — H{r)) in the other. For H{r) given in fll.3l) . we have 
(1 — H{r)) = (^_^_J-2q-) , while H{r) = -^jqi^^^; i.e., the same expression as (1 — H{r)) only 
with a in the latter replaced by —a. Therefore, if we consider the effective action with the 
background parameter a replaced by —a, the only change is that two potentials Vi^i+i/2{r) 
and V;+i/2,/(r) are interchanged, and so the two quantities ;_,_i^(oo) and S(^l_^_l^l^{oo) in 
f l3.23p are also interchanged. This shows that each partial wave contribution to the effective 
action with the background parameter —a is the same as the one with the parameter a 
and thus two effective actions with a and with —a have the same value. (This was true 
for the! classical action also). Similar behaviors, concerning the cases with a = ±1, were 
observed already in Ref . [l| . Based on this observation, consideration of the effective action 
for positive values of a is sufficient. 

The ODE system specified by ( 13.19p -( l3.2ip can easily be solved numerically. (With m = 
and a = ±1, analytic solution to this equation was found in In Fig. [2] we plot the 
solutions for a few cases. It clearly shows that the solutions approach constant values in the 
r ^ oo limit. In Fig. [3] we plot partial wave contributions with / = 0, 1/2, ■ ■ ■ for m = 1 
and a = 2. Note that P{1) ~ 0( j) when the angular momentum / becomes large. Since the 
degeneracy factor {21 + 1)(2/ + 2) is quadratic, this implies that rj<j^ in (13.221) behaves as 

in the large L limit. This divergent behavior will be canceled when we add the higher 
angular momentum contribution. 

D. WKB Calculation of the Higher Angular Momentum Part 

We now calculate the higher angular momentum part rj>j^, given in (12.61) . This can- 
not be computed numerically (as we have done for rj<j^) because very large partial wave 
contributions lead naively to a divergent result and require careful renormalization to en- 
sure a finite result. The large partial wave contribution depends on the regulating cutoff A, 
whose effect must be identified and isolated for renormalization; and this cannot easily be 
done numerically. However, this quantity (incorporating renormalization) can be calculated 
analytically in a WKB-type asymptotic series, assuming that the partial wave cutoff Jl is 
large enough. Here the higher angular momentum sum of the partial wave heat kernel. 
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S(i,i+i/2)ir) 



3.755 



3.750 



a=4 

Qf-3 



Q'=2 



a - 5 



3.745 



Qr=l 



(X 10^) 



3.735 



(a)Plot of the solutions 5';,;+i/2(f) for a = 1,2,3,4,5 (from the bottom) 
S(i+i/2,l){r) 



-3.765 
-3.770 
-3.775 
-3.780 
-3.785 



4 6 



(X 10') 



a-2 



a - 3 



a-4 



a-5 



(b)Plot of the solutions 5/+i/2.i(^) for a ~ 1,2, 3, 4, 5 (from the top) 

FIG. 2: Plots of Si^i^i/2{''') and 5'i+i/2,«('') when m = 1 and Z = 10. Note that the asymptotic 
values of Si^i/2,i{^) 'S'i,/+i/2(^) ^-^e roughly of the same magnitude but with opposite signs. 

J2j>.Jl ^ji^) with Fj{s) given by (12. 7p . may be described more explicitly by the form 

POO °° ^ _ _ 

Y: Fj{s) = j dr (2/+l)(2/ + 2){A(,,,+ .)(r,r;s)+A(,+ i,^.=,)(r,r;.) 

-A;5f(r,r;.)-A;;-.)(r,r;.)}. (3.26) 

Now, as explained in j^, we may use the y expansion for the modified radial proper-time 
Green function when / is large. When A(r, r; s) is expanded in terms with increasing number 
of derivatives of the potential, the scaling is such that when the / sum is approximated by 
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FIG. 3: Plots of partial wave contributions as a function of I for the cases with a = 1,2,3,4,5 
(from the bottom) and m = 1. 

the Euler-Maclaurin formula, this generates the large L expansion. For a generic radial 
potential V{r), this expansion has the following form: 



A(r, r; s) 



-sV{r) 



+ 



288^ 



1 + 

11 
360 



12 



6 



:s'V' 



^ ' 40^ ' 30 60 J 



10368 
43 



17{V')^V"s^ ^ s^iwyi' 



8640 



10080 



61 



15120 



{V 



+ 



2520 
1 



1008 



23 
5040 



(y(3))2^5 



19 
2520 



0(1 



(3.27) 



280 840 

Note that the terms are collected according to the total number of derivatives on V. [In 
addition to the terms already calculated in (3.16) of 2I, we have here included some higher 
order terms as well because they will be useful in finding the large Jl expansion of Fj>j^]. 
The large-L series expression for Tj^j^ can then be found by inserting (13.271) into fl3.26p 
with the generic potential V replaced by the potential Vij or V^^'^'^'^ in fl3.24p or fl3.25p . The 
summation over I in (13.261) can be performed using the Euler-Maclaurin summation method. 
The result is tantamount to the systematic WKB series, as was shown in Ref. j^. Since 
this procedure was described already in Appendix C of (2], we will not repeat it here. The 
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final result in the present potential can be presented as a series of the form 
Tj>j,= rfr|Q2(r)L2 + gi(r)L + Qiog(r)lnL + go(r) + Q-i(r)-4 



where 



Q2{r) 



%H{H -I) 



Qiog(r) 
Qo(r) - 



1 

4r 



6r(f2 + 4)7/2 



4ij2(iJ-l)2 + r2if'2j^ 

4(3f^ + 49f^ + 236f2 + 352)H^{H - if 



Q-iir 



-6(22P + 157P + 352P + 384)H{H - 1) + 16r(P + 5P + 4)(2if - 1)H' 
+r2(f2 + 4)2 {(3f2 + 8)i7'2 _ 4(2i7 _ 1)H"}] - Qi,^{r) In 
1 



^(f 2 _^ 4-) 9/2 



2(9P + 36P - 64P - 25Q)H\H - 1 



(3.28) 

(3.29) 
(3.30) 
(3.31) 



(3.32) 



i-6r^ + 25f^ + 368f^ + 128f^)H{H - 1) + 8rf2(f^ + 3r^ - A)(2H - l)H' 



(3.33) 



with f = The r-integration, with H{r) given in (11.31) . can be performed numerically. 
The next higher order terms, which are quite lengthy and so are not given here, can also be 
calculated in a straightforward manner. In fact we have calculated the quantity T jyj^ up 
to 0(-p)-terms, and the details of this calculation will be reported elsewhere 21 1. 



E. Results for the Total One-loop Effective Action 

Let us now put together the lower and higher angular momentum parts, f l3.22p and f l3.28p 
respectively, of the effective action, computed separately in the previous two subsections. 
Even if from the expression in (13.281) we keep only up to the terms of 0{LP) (i.e., up to the 
Qo-term in the integrand) and add it to the lower angular momentum part found numerically 
in the subsection C, we obtain a finite result in the limit of very large L, with divergent 
contributions from the two parts canceling each other. However, for moderately large L the 
resulting sum shows dependence on the cutoff value L. In other words, although the desired 
effective action should result if L is taken to be sufficiently large, the rate of convergence is 



quite slow (see Fig, 4(a)). This causes a problem in numerical efficiency. 
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FIG. 4: We have plotted L-dependence of the sum of lower and higher angular momentum parts. 
In (a) the upper (black) dots denote the case with all O ^-^^ terms ignored in Tjyj^, and slow 
convergence is evident, the lower (grey) squares the results after incorporating O ^-^j corrections, 
show better convergence. In (b) the (grey) squares, (purple) crosses, (blue) stars and (red) dots 
represent the cases obtained after we incorporate 8j^, -p- and corrections successively. 

One can accelerate the convergence by using the Richardson extrapolation method as 
described, say, in Ref. 1^, and this was in fact the method we used in the calculation 
of the instanton determinant in In this work we adopt a different, theoretically far 
more satisfying, approach to this problem. As we emphasized in SecJTll the effective action 
should not depend on the choice of our cutoff value L if the exact results for both 
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r(?7ip) 




mp 



FIG. 5: Plots of the modified eff'ective action as a function of mp. The (blue) dots, (black) stars, 
(brown) squares, (red) crosses and (purple) diamonds denote the values we get numerically for 
a = 1,2, 3, 4, 5 and the solid lines are for the associated large mass approximations. 

and rj>j^ were used. This implies that, leaving aside possible numerical inaccuracy in 
calculating rj<j^, the L-dependence in the sum for finite cutoff value L is really due to our 
ignoring —suppressed contributions in the WKB series (13.281) . So we can systematically 
improve the large L limit by adding higher order terms in the -^-series for the higher angular 
momentum part T j> . As we mentioned already, we have identified the -^-series up to terms 
of O(-p). When these higher-order terms of the -^-series are utilized, the situation changes 



dramatically: this is exhibited in Fig 4(b), From the figure we see that it is possible to 
achieve an L-independent result even for relatively lower and lower values of L if we include 
higher and higher order terms in for the large angular momentum part. Even L ~ 10 
produces good convergence. This reduces the number of numerical computations to be done 
in the low partial wave piece fl3.22p . Thus, use of the systematic WKB series for rj>j^ plays 
a pivotal role in reducing the computer time in our calculation. 

By the above procedure we have evaluated with high precision the renormalized effective 
action in given non-Abelian radial backgrounds. The results are shown in FigO In FiglJlwe 
plot (as a function of mp) our numerically accurate results for the modified effective action 
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r, defined in fl2.8p . with the choices a = 1,2,3,4,5 for the shape parameter in the radial 
profile function /(r) in (11.31) . We also plot, with solid lines, the predictions for the same 
quantity based on the large mass expansion (considered up to the order j;^^)- The large 
mass expansion is generally quite good when mp is large. But the validity range of the large 
mass expansion varies with the value of a; as a becomes large, the large mass expansion is 
reliable only when mp is significantly larger than the corresponding value of mp with small 
a. This is understandable, for the 'small' quantity in the large mass expansion is really the 
ratio of typical values of the fields and the derivatives of fields to the mass. The derivatives 
have larger values when a becomes large, and thus the ratio becomes small only when the 
mass has a comparatively larger value. 

TABLE I: Table of two parameters in (|3.34p for the five different values of a = 1, 2, 3, 4, 5. 
a 1 2 3 4 5 

Ai 0.145873 0.129759 0.113632 0.082787 0.037283 
A2 0.499(17) 0.305(81) 0.291(68) 0.292(11) 0.295(03) 



We now comment on the case with m = 0. In this case our numerical approach is not 
directly apphcable. Solutions to (I3.19p . S^i^i^i/2){r) and 5'(/+i/2,z)(r), behave like Inr and 
— Inr in the r ^ 00 limit, so the sum is finite. Analytic solutions to these equations, for the 
case of a = 1, can be found in jl|]. In this instanton case, the m limit was shown to be 
smooth, and furthermore to reproduce exactly 't Hooft's analytic massless result. We have 
found here that for other values of a, the m ^ limit is once again smooth, and from the 
numerical values of the effective action for m = jfo' ' ' ' ' iM' have found the following 
numerical extrapolations. In the case of a = 1 the extrapolation is 

r(m) = + ^2 In m 

= 0.145873(29) + 0.499(17)m2 In m (3.34) 

In this case, analytic expressions for two leading coefficients in the small mass expansions are 
:mown, with the resuhs: Ai = + i(l - In 2) - 2C'(-1) ~ 0.14587331 ■ ■ ■ Q; and ^2 = ^ 
9], [22I. Note the remarkable agreement between these analytic results and the numerical 
extrapolation in ( 13.34p . Numerically determined (through extrapolation) values of Ai and 
A2 for the cases with other values of a are presented in Table [H We have not yet found 
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simple analytic expressions for these leading small mass terms, although we suspect this 
should be possible. 



IV. QUASI ABELIAN BACKGROUND FIELDS 

In this section we consider the cases of quasi Abelian background fields, where an Abelian 
gauge field is embedded in the SU(2) Yang-Mills gauge fields. Using an appropriate isospin 
rotation we can take the non-vanishing components only in the 3rd direction. Then the 
vector potential can be written as in (11.21) . and we take the radial profile function g{r) in the 
step-like form (II. 4p . Recall that there are three arbitrary parameters /?, £o; and B. When 
P = it corresponds to the case of uniform field strength studied in [2]. When P ^ 0, 
this potential describes a bubble shape with a radius Rq = ^q/^/B, and the dimensionless 
parameter (3 is related to the thickness of the bubble. Illustrative graphs of the radial 
function g{r) for various values of Rq are drawn in FiglHl 

The field strength tensor of the quasi Abelian gauge fields is 

F^^ = -2r]^^iu'g{r)T^ + —{x^rj^xi - x^r]^xi)u'g\r)T'^ , (4.1) 
The classical action can be evaluated as 



]^jct^xtiFl^ = Att^ drr'^[8g{rf + 4rg'{r)g{r)+r^g'{ry^ (4.2) 

(PolyLog(3,-e2^«)-PolyLog(5,-e2«")). (4.3) 
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In the large Rq limit, this quantity can be approximated as 

207r2 /4 10 , ^ „, 1 



, (fxtiF' ^ ( -R^r. + — (tt^ - Q)Rf. + — (Ttt^ - 60'k^)Ro + . . . . (4.4) 
2 J '^'^ V3 ° 9 ^ ^ ° 36^ ' J ^ ' 

We turn to the evaluation of the one-loop effective action. When the background function 
is given in the form (11.40 . based on dimensional considerations, we may cast the renormal- 
ized action into the form (12. 9p . and so simply compute T(^m/\/B^. Furthermore, since 
f {m/^/B^ is a function just of we may set -B = 1 hereafter, and it can be restored 
later. 
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(b)Plots of g{r) with fixed i?o = 3 



FIG. 6: Plots of the function g{r) (in units of B) for various values of i?o and /?. In (a), plots are 
drawn for Rq = 1, 2, 3, 5, 10 with fixed (3=1. In (b), plots are drawn for (3 = 5,2, 1, 1/2, 1/5 when 
i?o = 3 

A. Large Mass and Derivative Expansions 

In the large mass expansion as given by fl2.1ip . the modified effective action r(m) can be 
approximated by 

f Mlm = fffilnm + fSi^ + f^i^ + ■ ■ • , (4.5) 

where 

f a = -^j^ dr^' + '^rg{r)g'{r) + rV(r)') (4.6) 
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720 Jo 



.(4) ' 



10080 7o 



2Ar^g{r) (l5g'{r) + r (9/(r) + rg^^\r) 
(221g'{rf + 9r^g"{rf + 2r^'(r) (7l/(r) + 6r^(^^(r)))] (4.7) 
- -2688r^^(r)^ - 540^(r)^'(r) - 2688r^^(r) V(^) + 595rg'{rf 
-U56r^g{rfg'{rf - 392r^g{r)g'{rf - 49rV(^)^ + 540rg{r)g"{r) 
+U80r^g'{r)g"{r) + 1837rV'(^)^ + 1620r^ g{r)g^^\r) 
+2356r^ g'{r)g^^\r) + 830rV(^)^^^n^) + S7r^g'''^\rf + 504r^g{r)g^^\r) 
+380r^g'{r)g^'^\r) + 52r^ g" (r) g^'^\r) + 18r^ (2^(r) + rg'{r)) g^^\r)] (4.8) 

For g{r) given by the form (ll.4p it is straightforward to evaluate numerically the integrals 
in dUD, (gTD and fO|) . 

There is another well-known approximation method - the derivative expansion . The 
leading term in this expansion is given in (12.141) . where ±iEi and ±iE2 are four different 
eigenvalues of the matrix F = {F^^), and in this quasi- Abelian case, can be expressed in 
terms of the function g{r) as 

E, = 2g{r) , E, = 2g{r) + . (4.9) 

dr 

For the derivative expansion, one may thus insert the results (14.91) into our formula (I2.14p and 
evaluate numerically the r and s integrals. This is also straightforward. These large mass 
and derivative expansion approximations will be compared with the exact results below. 



B. Numerically Accurate Calculation of the Lower Angular Momentum Part 

When the gauge vector fields have the form (11.21) . the partial waves are specified by the 
quantum numbers J = (/, Z3, ts, Ts). The radial Hamiltonian becomes [2I 

n^„t, ^ -^(M3,t3) = -^(0 + Hr)ht3 + r^9{rf (4.10) 

Note that this Hamiltonian does not change when we simultaneously change the sign of 
and l^. Using this symmetry we can set ^3 = 1/2 without loss of generality. Then the first 
(lower angular momentum) part of the renormalized action can be written as 
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11. 



FIG. 7: Plot of solutions to (fiJ3]) with /a = 4, 3, • • • , -3, -4 (from the top) when ^ = 4, /3 = 1, 
Rq = 3, and m = 1. Note that the curves are essentially flat beyond r > Rq. 

As in the case of Sec JIIII the ratio of two determinants is determined by the asymptotic value 
of the function 



Sr 



1 i^i,hir) 

which satisfies the following differential equation 



'dS, 



dr"^ \ dr 
with boundary conditions 



dr 



+ {- + 2m:^^^±i^ > "^^^'^''^ 
\r l2i+i{mr) 



%3)(0) = 5(M3)(0) = 0. 



The potential term ?/(/,/3)(?") in (14.131) is given by 



(4.12) 



(4.13) 



(4.14) 



U^i,i,){r)=4kg{r)+r\g{r)f 



(4.15) 



This differential equation (14.131) can be solved numerically. For the case with f3 = 1 and 

, 3, 4 when / = 4 in FigJTl 
/, we get the partial wave 



-Ro = 3, we plot the solutions of (14.131) with Z3 = —4, —3, 
For a given /, after summing all the values with = —I,- 
contributions: 

^(0 = E ^(/,/3)(oo). 

For large / we find P{1) behaving like j. 



(4.16) 
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C. WKB Calculation of the Higher Angular Momentum Part 



To compute the renormalized effective action we must identify the higher angular mo- 
mentum part. The leading contributions in the large angular momentum limit is written in 
the form (13.281) . We present some of them in explicit forms: 

(4.17) 



Q2{r) 
Qi{r) ^ 
Qiog{r] 



2r3(3f2 + 8) 2 

(f2_^ 4)3/2 3 5 
3 



(4.18) 
(4.19) 



Qo{r 



90(f2+ 4)7/2 



12r^(5f^ + 28r^ + 32)^^ - 30 (9f^ + 47f^ + 40f^ + 64)^' 

+20r(3f^ + 32f^ + 88^^ + 32)^^' + 5r^{f^ + Af {{3f^ + 8)^'^ - 8^/}" 
fir \ 



-Qiogir) In 



g^4~4(^~2 + 4)^4 _ (^4~8 ^ 7~6 ^ 4g~4 ^ 1x52^2 ^ 1024)^' 



(4.20) 



Q-i{r 



with r = Higher order terms will be presented in [2]J. 



4(f2 + 4)9/2 

16r(f6 + 15f^ + 52f2 + 32)^^' - ir^i^r^ + 4f {{r^ + 2)g'^ + f^^/}] , (4.21) 



D. Results for the Total One- loop Effective Action 

As in the non-Abelian case treated in Sec lIIIt the large L divergence of the numerical 
results for the low partial wave contribution is canceled by the large L divergence in (I3.28p . 
found analytically from radial WKB. The combination of these lower and higher angular 
momentum parts is then L- independent in the large L limit, as we have seen in the last 
section. In Fig. [8|, we plot the effective action for various values of the parameters m, Rq 
with fixed value of (3 = 1. 

In Fig. [9l we plot the effective action for various values of the parameters m, (3 with fixed 
value of -Rq = 0. Its large mass expansion and its derivative expansion are drawn together for 
comparison. Clearly, the derivative expansion becomes more accurate when the parameter 
/3, which represents of the derivative scale of the background function g{r), becomes smaller. 
Note that the result of leading derivative expansion with i?o = becomes independent of /3 
when it is multiplied by Further note that the derivative expansion is generally much 
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(c) (d) 

FIG. 8: Plots of the effective action for different Rq values, i.e., i?o = in (a), i?o = 3 in (b), 
i?o = 5 in (c), and Rq = 10 in (d), assuming (3 = 1. The solid line in each figure denotes the result 
of the derivative expansion of the effective action while the dashed line denotes the result based 
on the large mass expansion. 

better than the large mass expansion for smaller values of m. This is despite the fact that 
we have just used the very leading order of the derivative expansion fl2.14p . in which we 
simply take the Euler-Heisenberg constant field result, and then replaced the constant fields 
by their inhomogeneous forms in the effective Lagrangian. The superiority of the derivative 
expansion is because the derivative expansion is in fact a resummed version of the large 
mass expansion - the large mass expansion is an expansion both in powers of the field and 
in derivatives of the field, while the derivative expansion is just an expansion in derivatives 
of the field, with all terms in the large mass expansion not involving derivatives having been 
resummed. This difference is clearly reflected in the plots shown in Fig. [HI 
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FIG. 9: Plot of the effective action multiplied by (3^ as a function of m when Rq = 0. The (blue) 
dots, (red) diamonds, (grey) squares and (brown) crosses denote the numerical values for /? = |, 
P = ^, (3 = 1 and [3 = 2, respectively. Each of dashed lines denotes the associated large mass 
expansion. The solid line denotes the leading derivative expansion which is independent of (3 when 
i?o = 0. 

V. CONCLUSIONS 



To conclude, we have presented explicit computations of the renormalized one loop ef- 
fective action for gauge field backgrounds that possess a radial symmetry, such that the 
associated spectral problem can be decomposed into partial waves. We considered one class 
of non-Abelian backgrounds, and another class of quasi- Abelian backgrounds, each charac- 
terized by a radial profile function. The computation has been performed using the partial 
wave cutoff method developed in [l, 2], and we have further refined the numerical efficiency 
and precision by including higher order terms in the analytic radial WKB expression for 
the large partial wave contribution. The main conclusion is that the method works very 
efficiently and simply. With the incorporation of these higher order analytic terms, the nu- 
merical part of the computation is simplified because we do not need to take such high partial 
waves in the numerical computations. The method is not much more complicated to imple- 
ment than the derivative expansion or the large mass expansion, and is much more accurate, 
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especially in probing the small mass region. So, we can now reliably compute the renor- 
malized determinant of any fluctuation problem whenever there is a radial symmetry. This 
method has now been tested successfully in gauge theories [l, 2| and in self-interacting scalar 



fleld theories 
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26| . The physical renormalization conditions are quite different in 



these various theories, but the WKB analysis correctly encodes the renormalization physics 
in each case. The method has also motivated a recent extension of the Gel'fand-Yaglom the- 



orem (for the determinant o 



ordinary differential operators), to partial differential operators 



that are radially separable 2?! . 

An important generalization is to fermion flelds. This can be done by converting the 
fermion problem to second order form, and then using the scalar method we have described 
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30, 



here. Aspects of this idea have been addressed in various approximation schemes 
311 ]: however, it should be possible to develop a more direct and numerically exact fermionic 
approach, along the lines of the scalar partial wave cutoff method described here. Finally, 
while the method is restricted to backgrounds for which the fluctuation spectral problem is 
separable into partial waves, this includes a relatively large class of physically interesting 
cases, such as vortices, monopoles, and instantons. 
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